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Abstract 

We study the astrophysical reaction rate for the formation of ^Be through the three 
body reaction a{an,'y). This reaction is one of the key reactions which could bridge 
the mass gap at A = 8 nuclear systems to produce intermediate-to-heavy mass elements 
in alpha- and neutron-rich environments such as r-process nucleosynthesis in supernova 
explosions, s-proccss nucleosynthesis in asymptotic giant branch (AGB) stars, and pri- 
mordial nucleosynthesis in baryon inhomogeneous cosmological models. To calculate the 
thermonuclear reaction rate in a wide range of temperatures, we numerically integrate the 
thermal average of cross sections assuming a two-steps formation through a metastable 
^Be, a 4- a v=^ ^Be(n,7)^Be. Off-resonant and on-resonant contributions from the ground 
state in *Bc arc taken into account. As input cross section, we adopt the latest ex- 
perimental data by photodisintegration of ^Be with laser-electron photon beams, which 
covers all relevant resonances in ^Be. Experimental data near the neutron threshold are 
added with 7-ray flux corrections and a new least-squares analysis is made to deduce 
resonance parameters in the Breit-Wigner formulation. Based on the photodisintegration 
cross section, we provide the reaction rate for a{an, 7)^Be in the temperature range from 
T9=10~^ to T9=10^ (Tg is the temperature in units of 10^ K) both in the tabular form 
and in the analytical form for potential usage in nuclear reaction network calculations. 
The calculated reaction rate is compared with the reaction rates of the CF88 and the 
NACRE compilations. The CF88 rate, which is based on the photoneutron cross section 
for the 1/2+ state in ^Be by Berman et al, is valid at Tg > 0.028 due to lack of the 
off-resonant contribution. The CF88 rate differs from the present rate by a factor of two 
in a temperature range Tg > 0.1. The NACRE rate, which adopted different sources of 
experimental information on resonance states in ^Be, is 4-12 times larger than the present 
rate at Tg < 0.028, but is consistent with the present rate to within ±20% at Tg > 0.1. 
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1 Introduction 



The interplay between nuclear physics and astrophysics is essential to clarify the origin of 
elements in the Universe. Scenarios of the nucleosynthesis in stars rely on the determination of 
important nuclear reactions and structures along the reaction path A long way to create 
heavy elements starting from light elements (p, n, a) has to pass through the gaps in mass 
number at A=5 and A=8. The bridge to heavy elements such as ^^C invokes slow three-body 
reactions and then links to heavier elements up to iron elements and beyond. The triple alpha 
reaction for the formation of ^^C is usually the key reaction in stars at the helium burning 
stage. 

In a neutron-rich environment, however, there is an important three-body reaction. If 
neutrons are abundant enough together with alphas, the formation of ^Be can proceed through 
the reaction a{an, 7). This reaction followed by ^Be(a,n)^^C can compete with the triple alpha 
reaction and may dominate over, depending on the astrophysical environment. 

The importance of the aan reaction has been pointed out in the r-process in the hot 
bubble of a core-collapse supernova 0, ^. After the core bounce in the supernova, a nascent 
neutron star is formed at the center. The outer material of this proto-neutron star is heated 
by intense neutrinos emitted from the proto-neutron star's inside and is ejected as a wind. 
A high-entropy bubble in this neutrino-driven wind is thought to be an ideal site for the r- 
process nucleosynthesis. Because of the high temperature and the low density, the matter is 
dissociated into neutrons and protons once, and they reassemble into light elements, mostly 
alphas. The nucleosynthesis starts from alphas and remaining neutrons when the temperature 
decreases below Tg ~ 5 in the expanding matter. Since the triple alpha reaction is too slow 
at the temperature-density conditions in the hot bubble, the aan reaction forms a bypass 
flow connecting to ^Be(a,n)^^C. Therefore, the determination of the reaction rate for the ^Be 
formation is essential for the initial condition of the r-process, namely, the amounts of seed 
elements and remaining neutrons. Because of the high entropy (i.e. low density), the three 
body reaction proceeds slowly, resulting in seed elements with a high neutron-to-seed ratio for 
a successful r-process up to A=200. 

In recent studies of the r-process in the neutrino-driven wind, the importance of this reaction 
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is more evident. It was shown that the r-process can be successful in the wind flow with a very 
short expansion time and a moderate entropy Q . In this rapid expansion, there is essentially 
no time for the triple alpha reaction to set in and, instead, the aan reaction becomes again 
crucial for the r-process. The necessary expansion time is determined in comparison with the 
reaction time scale of the aan reaction [Q. A more recent study reveals that the reaction 
flow can go through the neutron-rich region of light elements in the rapidly expanding wind 
1^, 1^. Therefore, the competition of the aan reaction with other bypass reactions in the short 
expansion time should be carefully considered. 

The aan reaction may play an important role in other astrophysical environments. In the 
systematic study of alpha-rich freeze-out, the aan reaction may regulate the production of 
elements up to A=100 even if it does not lead to the r-process [§,0. In the neutron-rich He 
layer of low-mass asymptotic giant branch (AGB) stars, the aan reaction might be able to take 
place below T9=0.3. This may change the amount of neutrons to determine the abundance 
of s-process elements [|lOl. In the inhomogeneous Big-Bang nucleosynthesis, the aan reaction 
might contribute to the production of heavier elements beyond A=8 in neutron-rich zones at 
temperature around T9=0.1 ||rT|. 

Most of nucleosynthesis studies treating the aan reaction adopt the astrophysical reaction 
rate from the compilation of CF88 |T^. A more recent evaluation of the ^Be formation rate 
can be found in the NACRE compilation of the thermonuclear reaction rates . 



Since the direct measurement of the a(an, 7)^Be reaction is impossible, the experimental 
efforts to provide information on excited states in ^Be have been made (see [|l3l , Q and the 
references therein). A ternary process hardly plays a role in the formation of ^Be. Instead, the 
reaction reflects nuclear structure of ^Be with large neutron-decay widths of low-lying states in 
^Be. The reaction is described by a + a ^ '^Be(n,7)^Be, which takes place during the lifetime of 
^Be (~ 10~^^ s) in such suitable astrophysical site as the neutrino-driven wind. Therefore, the 
resonance structure of ^Be near the threshold for '^Be+n at 1.665 MeV is of great importance 



| 15| . It is noted that the reaction a + n ^ He(a,7) Be with the lifetime of He (~ 10 s) 
can make an additional contribution to the formation of ^Be |1T6|, [17 . 



Recently, a systematic measurement of the Be(7,n) Be reaction was made with laser- 



4 



electron photon beams in the energy range of astrophysical relevance |]T^. This study has 
provided photoneutron cross sections as a function of 7-ray energy. The data enables us to 
derive the neutron capture cross section of ®Be with the detailed-balance theorem, which is 
necessary in the calculation of the reaction rate for the ^Be formation. This systematic treat- 
ment of all the resonances is preferable in contrast with the preceding treatment which adopted 
different sources of experimental information on the resonances together with some assump- 
tions. 

The brief report |I^ provided the reaction rate of < aan > as a function of temperature. 
However, the reaction rate was not tabulated, but shown only graphically as a ratio with 
respect to the NACRE compilation. The temperature range was limited between T9=0.1 and 
T9=10, where only the on-resonant contribution from ^Be was taken into account. In contrast, 
the NACRE compilation treated a broader temperature range down to T9=10^'^ for possible 
astrophysical applications. At temperatures below T9=0.028, the effective energy window for 
the aan reaction is far below the resonance energy of ®Be. In this case, the off-resonant 
contribution for the formation of ^Be from a + a becomes important. This contribution was 
first shown to be important in the determination of the triple alpha reaction in accreting white 



dwarfs and neutron stars 18 



In this paper, we calculate the reaction rate for the ^Be formation taking into account both 
on- and off-resonant contributions based on the latest experimental information on photoneu- 
tron cross sections for ^Be. It is to be noted that we add data points immediately above the 
neutron threshold since the brief report |0. We provide the reaction rate in the wide range 
of temperatures from T9=10~'^ to T9=10^ in the analytical fitting formula as well as in the 
tabular form. 

We arrange this article as follows. In section ^ we briefiy review the preceding efforts to 
derive the aan reaction rate from the point of view of experimental information and exper- 
imental methods adopted. In section |^, we describe the photoneutron cross section for ^Be 
adopted to derive the reaction rate for the neutron capture of ®Be. The description of the 
experimental method is limited to minimum unless it is associated with added data and a new 
analysis. In section 11, we explain the theoretical method to derive the aan reaction rate from 
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the experimental data. In section ^, we provide the numerical results of the reaction rate in a 
tabular form and plots. We also provide the fitting formulae of the reaction rate. We compare 
the result with the reaction rate in the CF88 and NACRE compilations. The article is closed 
with a summary in section |^. 



2 Resonances of ^Be 

Fowler, Caughlan, and Zimmerman [|l^ presented the reaction rate of a (an, 7)^Be in their 
second volume of Thermonuclear Reaction Rates published in 1975 (hereafter referred to as 
FCZ II). The rate has not been revised in the third volume (1983) of the series and 
subsequently passed on to the two compilations of Caughlan, Fowler, Harris, and Zimmerman 
(1985) [ig and of Caughlan and Fowler (1988) ^ (hereafter referred to as CF88) without 
modification. It was nearly a quarter century after the FCZ II that Angulo et al. |T^ updated 
the reaction rate (hereafter referred to as NACRE). 

Early measurements of photoneutron cross sections for ^Be were carried out with monochro- 
matic 7 rays from radioactive isotopes p2| , ^ [2^ , p5| , ^ and with bremsstrahlung ^ 



The data of Herman et al. exhibit the maximum cross section 1.6 mb at 6 keV above the 
n + ^He threshold (E^ = 1.665 MeV) and decrease to 1.2 mb at 40 keV ||2^. Parametrizing 



the photoneutron cross sections by cr^„ = 1.6 exp(-4E) mb (E in MeV) with E = E^ - Ef [jT3 



resulted in the CF88 reaction rate N^< aan > as 



Nl {aan) = 2.59 x l^-'^^^^^^^^Wv,, (1) 
' [T|(l + O.344T9)] ' ^ ' 

where Tg gives temperature as T = TgX 10^ K and N^ is Avogadro's number. The low energy 
experimental data were recently investigated in a semi-classical model to derive the a(an, 7)^He 
reaction rate |]29[ . 



Woosley and Hoffman discussed the reaction rate in a narrow level approximation for two 
lowest states (the 1/2+ and 5/2^ states) in the n + ^Be channel. It was pointed out that when 



the resonance parameters from the (e,e') data of Clerc et al. are employed (i.e., E^i = 1.78 
± 0.03 MeV and F^i = 0.30 ± 0.12 eV for 1/2+ and E^2 = 2.44 ± 0.02 MeV and F^s = 0.089 
± 0.01 eV for 5/2^), the rate is consistent with the CF88 rate to within 30 % in the range Tg 



= 1 - 3 [|]. The narrow level approximation holds for the 5/2 state with F = 0.77 ± 0.15 keV, 
but not for the 1/2+ state with F = 217 ± 10 keV [|15]. In the discussion, the discrepancy in 



Fjxi was also pointed out between Ref. and the latest compilation of nuclear physics data 
| 15| which adopted E^^i = 1.684 MeV ± 7 keV from a different (e,e') data of Kuechler et al. 

After the FCZ II, some new measurements of a^n were carried out with both monochromatic 
7 rays from radioactive isotopes and with bremsstrahlung [^, |3^. An R- matrix fit to the 
a^n of Fujishiro et al. |32[ in one level approximation gave parameters B = LOGlg^g mb, eji 

-+0.019 o2fwi2 



192 keV, and = 67 keV [gg which can be converted to B(E1 i) = O.lOGj:^;^}^ e^fm^, E^ 
1.733 MeV, and F = F„ = 227 keV, respectively. The last quantity was evaluated from 



F„ = 2Jeji{Ej — Et) with — Et = 67 keV. A major difference between the photonuclear 



reaction and the electron scattering is that the former B(E1 1) |32| is a factor of two lager than 
the latter one (B(E1 i) = 0.050 ± 0.020 |30| and B(E1 i) = 0.054 ± 0.004 [P). According to 
the Barker's analysis, the main reason for the discrepancy may be due to a potential problem 
of background subtraction in the electron scattering. Barker pointed out that the (e,e') result 
may represent not a full but a partial strength of B(E1 J.). 

The NACRE rate of Angulo et al. is based on three states in ^Be, i.e., the 1/2+ , l/2~, and 
5/2+ states, neglecting the 5/2^ state. Resonance parameters (E^,, F„, and F^) employed for 
the 1/2+ state were taken from the cr^„ measurements of Fujishiro et al. ||3^, where F^ was 
0.51 ± 0.10 eV. For the broad 1/2" state at E^ = 2.78 ± 0.12 MeV with F = 1080 ± 110 keV 
||T5| , F^ = 1 W.u. (0.45 eV) for an Ml decay was assumed with a 80 % uncertainty [^. For 



the 5/2+ state at E^. = 3.049 MeV ± 9 keV with F = 282 ± 11 keV pl, F^ = 0.90 ± 0.45 eV 
which is much larger than the (e,e') result (0.30 ± 0.25 eV) of Clerc et al. was adopted with 
the increased uncertainty based on the ct^„ data of Hughes et al. PB| . 

Recently, Utsunomiya et al. made a new measurement of a^n for ^Be in the range of E^ = 
1.78 - 6.11 with quasi-monochromatic 7 rays produced by means of inverse Compton scattering 
of laser photons [|T4|. The data essentially include all states of astrophysical interest. A least- 
squares fit to the data was carried out to extract the Breit-Wigner resonance parameters. 
Using the (7,n) cross sections, the a (an, 7)^Be reaction rate was evaluated. This evaluation 
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was based only on the resonant contribution to the reaction a + a ^ ^Be with = 6.8 ± 1.7 
eV. Consequently, the reaction rate is valid only at Tg > 0.028 as is the case for the triple alpha 
reaction |jl8[ . 



The 5/2^ state was treated differently in the NACRE rate and the rate of Utsunomiya et 
al. It was completely omitted in the former, whereas it was taken into account in the latter 
assuming its 100 % decay to the '^Be(O) + n channel. The 5/2^ state in ^Be was investigated in 
the '^Li(^He,p)^Be reaction |^ and in (3 decay of ^Li and found to decay 7.5 ± 1.5 % and 
6.4 ± 1.2 % to ^Be(O) + n, respectively. It is reported that the 5/2^ state decays dominantly 
to 5He(0) + a |9[. The 5/2" state with V = 0.77 ± 0.15 keV lies at = 2.4294 MeV ± 1.3 
keV below the ^He + ^He threshold energy (2.467 MeV). However, since the ground state in 
^He has a large width (0.60 ± 0.02 MeV) the state can decay into the low-energy tail of 
the ground state. 

A potential contribution of the ^He(O) + a channel to the reaction rate was investigated 
by Buchmann |jl6| where a 5 % branching of the 5/2~ state to the n + ®Be decay channel was 
taken. Summing up all possible contributions through the ^He + a channel not only from the 
5/2~ state but also from the 1/2+ and 1/2^ states resulted in a reaction rate that is one order 
of magnitude smaller than through the n + ^Be channel at high temperatures, Tg > 2 (Fig. 9 
of Ref. |16l). 

In a recent experiment of /^-delayed particle decay of ^C, a large background continuum 



in B which decays to the Li channel was observed above 2.9 MeV [Q. An R-matrix fit 
revealed that a low energy part of the continuum consists of a tail of the ground state 0]. 
The charge-symmetry analogue of this background continuum in ^Be may contribute to the 
a{an, 7)^Be reaction rate. An attempt to estimate the contribution of the analogue background 
continuum to the a{an, 7)^Be reaction rate through the ^He + a channel was made, showing 
that it possibly dominates the rate for high temperatures above Tg ^ 4 |1^, This result 
may be qualitatively consistent with the result of a recent microscopic three-cluster study of 
photoneutron cross sections for ^Be which showed that the ^He + a channel is of growing 



importance above E^ ^ 4 MeV [^ . 

In the present paper, we evaluate the a (an, 7)^Be reaction rate through the ®Be -|- n channel. 
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including a fractional contribution of the 5/2^ state. The contribution through the ^He + a 
channel is not included because of a lack of sufficient experimental information on the ^He + 
a channel. Clearly, the evaluation of the stellar reaction rate at the next level must await a 
new experimental outcome for the decay of relevant states in ^Be through the ^He + a channel 
(particle- and 7- widths and branching ratios of decays) which is expected in the proposed 
experiment of /3-delayed particle decay of ^Li [|16], |17| . 



3 Photodisintegration of ^Be 
3.1 Experiment 

We limit a description of the experiment to minimum and concentrate on new facets of the data 



and data analysis. Many of the details can be found in |]I4|. The experiment was performed 
at the National Institute of Advanced Industrial Science and Technology (AIST) which is 
reorganized from the Electrotechnical Laboratory. Quasi-monochromatic 7-ray beams were 
generated by means of inverse Compton scattering of YLF: Nd laser photons with A = 1053 
nm (1.2 eV) from relativistic electrons stored in an accumulator ring TERAS. The energy of 
the laser-electron photon beams (LEPBs) was varied in the energy range of 1.69 - 6.11 MeV 
by changing the electron beam energy from 308 to 587 MeV. The LEPBs were collimated into 
a spot 2 mm in diameter with a 20 cm thick Pb block at 550 cm from the head-on collision 
point. The beam intensity which gradually decreased with the electron beam current was in 
the range of 10^ - 10^ cps. A 4 cm thick ^Be was irradiated and neutrons were detected with 
four BF3 counters embedded in a 30 cm polyethylene cube. A BGO detector was used as a flux 
monitor of the LEPB. 

Experimentally two reaction channels of n + ^Be and ^He + a cannot be distinguished. 
Different cross sections may result depending on the reaction channels because the neutron 
detection efficiency may be different in the two channels. This complication could occur for 
the 5/2~ state which decays to the two channels, whereas other states are free from such 
complication because they are known to predominantly decay to the n + '^Be channel [0. The 



average energy of neutrons emerging from the decay of the 5/2 state to the two channels was 
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calculated. The energy of this state is well defined at 2.4294 MeV (± 1.3 keV) with T = 0.77 
keV (± 0.15 keV). When this state decays to the n + '^Be channel at 1.665 MeV, the average 
laboratory energy of neutrons is 382 keV. On the other hand, this state can decay to the ^He 
+ a channel at 2.467 MeV with the low-energy tail of T (^He) = 0.60 MeV (± 0.02 MeV). The 
average energy of neutrons resulting from ^He was calculated to be 156 keV. 

Figure |l] shows the detection efficiency of the neutron detector as a function of neutron 
energy. The efficiency was calibrated with a 265 keV-neutron source of ^^NaOH + D2O and a 
^^^Cf standard source. The energy dependence of the efficiency was calculated with a Monte 
Carlo code MCNP |3[. The efficiency was 6.62 % at 382 keV and 6.65 % at 156 keV. Thus, the 
difference in detection efficiency between the two reaction channels is negligible. Cross sections 
presented in this paper are based on the detection efficiency for the n + ^Be reaction channel. 

Figure |^ shows an energy spectrum of the LEPB measured with a high-purity Ge detector. 
The LEPB has a characteristic low-energy tail that is determined by the electron beam emit- 
tance, the laser optics and the collimator size. When the low-energy tail crosses the neutron 
threshold of ®Be + n (1.665 MeV), not the total but a partial flux of the LEPB is responsible 
for photodisintegration of ^Be. This occurs for LEPBs whose peak energies are very close to the 
neutron threshold energy. Previously photoneutron cross sections were presented at energies 
above 1.78 MeV without flux corrections. As a result, cross sections in the peak region of the 



1/2+ state were missing in [|T4 



The flux correction was made for LEPBs by means of a Monte Carlo simulation. Energy 
spectra of the LEPBs taken with the Ge detector were reproduced with a Monte Carlo code 
EGS4 ||4^ by simulating electromagnetic interactions between a LEPB and the Ge detector. 
The energy distribution of the LEPB which best reproduced the response of the Ge detector is 
shown by the solid line in Fig. ^. 

Figure ^ shows photoneutron cross sections obtained in the present experiment. The flux 
correction was made for the data in the energy range from the neutron threshold to 2 MeV. 
These data are plotted at the average energies by the solid circles in the flgure. Note that the 
three data points at 1.69, 1.71, and 1.73 MeV which covered the peak region of the 1/2+ state 
are consistent with the data of Fujishiro et al. It was essential to have these data points in a 
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new analysis because the statistical uncertainty of the Fujishiro data is quite large. 

Including all the data with the flux correction, a new least-squares fit was performed to 
extract the Breit-Wigner parameters for resonance states in ^Be. The same fitting procedure 
was followed using a straight-line cross section (cr = 0.38 - 1.21 [mb], in MeV) 



as m 



for higher-lying states in the energy region up to 6.1 MeV. The El parametrization was used 
for the positive-parity states (1/2+ and 5/2+), while the Ml parametrization for the negative- 
parity state (5/2~). Due to the fact that the energy spread of the LEPBs was much larger than 
the total width (< 1 keV) of the 5/2^ state, the line shape of the state was not determined 
experimentally. As a result, only was deduced from the integrated cross section. The best 
fit is shown by the solid lines. The best fit parameters are listed in Table 1. 

It is to be pointed out that there are large differences in B(E1 1) and B(M1 1) for ^Be 
states between photodisintegration and electron scattering. While the present B(E1 1) for the 
1/2+ state is consistent with the result of the photodisintegration experiment with radioactive 
isotopes (0.106l|];[ji^) |]32[| , it is larger by a factor of two than the results of electron scattering 
(0.050 ± 0.020 g and 0.054 ± 0.004 [§T)). Similarly, the present B(Eli) for the 5/2+ state 
is about 4 times larger than that (0.010 ± 0.008) from (e, e') [^, whereas the B(M1 1) for the 
5/2- state is about half that (0.536 ± 0.060) from (e, e') [0. 

Regarding the difference for 1/2+, Barker previously pointed out a potential problem of 
background subtraction in electron scattering. As can be seen in the line-shape (Fig. |^), the 
1/2+ and 5/2+ states overlap each other with substantial tails underneath the narrow 5/2^ 
state. The electron scattering strongly enhanced the excitation of the Ml state (5/2~) adjacent 
to the El states (1/2+ and 5/2+) 0, p. The authors of Ref. |g reahzed a background 
subtraction problem in the discussion of the 1/2+ state (Section 5.3 of Ref. [30|), with the 
statement "Because of the large background and the limited statistical accuracy of this small 
peak, the contribution of the long tail to be expected from the (7,n) cross section is not 
included in the result. Therefore, it is reasonable to make a comparison with the (7,n) cross 
section integrated up to 2 MeV only." Indeed, the B(E1 1) resulting from the present (7,n) 
cross section integrated up to 2 MeV is in good agreement with that given in the electron 
scattering. Similarly, the B(E1 1) for the weakly populated 5/2+ state can be underestimated 
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in the electron scattering for the background problem. This explains, however, the difference 
only qualitatively which is larger for the 5/2+ state (4 times) than for the 1/2+ state (2 times). 

Regarding the differences for the 5/2^ state and the 5/2+ state, a broad Ml state (l/2~), 
which supposedly exists over the 2.5 - 3.5 MeV region |3^, complicates the situation. The 
present least-squares fit did not single out the l/2~ state which was reported to be a broad 
state at 2.78 MeV with F ~ 1 MeV [|l^]. If one treats the l/2~ state exactly in the same manner 
as in the NACRE compilation (E^ = 2.78 ± 0.12 MeV, T = 1080 ± 110 keV, = 1 W.u. 
(0.45 eV)), the peak cross section at 2.78 MeV is found to be 0.066 mb. This cross section is 
one order of magnitude smaller than the experimental cross section in the corresponding energy 
region. It is interesting to note that the present for 5/2+ (1.24 ± 0.02 eV) is accidentally 
close to the sum of the F^'s employed in the NACRE compilation for 5/2+ (0.90 ± 0.45 eV) 
and 1/2^ (0.45 eV) to within the associated uncertainty. A more quantitative understanding 
of this energy region is necessary. 

3.2 Photoneutron and neutron capture cross sections 

We here present analytic formulae to express the cross sections with the resonance parameters 
which are determined by the least-squares fit described in the previous section. 
The photoneutron cross sections for each resonance state of ^Be are given as, 

, , y - [181. 9(^^ - 1665.3)1^/2 

cj^n 1/2+ = 6.654 X IQ-^E^j-^ ^ mb , (2) 

^'^ ' ' ^[(E^- 1735.1)2 + 181.9(^^- 1665.3)]^ ^' 

274 72 

a^„(5/2+) = 7.805 x IQ-^E^-^ — [mb], (3) 

^'^ ' ' ^(E^- 3077.0)2 + 75471.^ J' v J 

a^J5/2-) = 3.4616 x lO^E^"^-^ \- [mbl. (4) 

^'^ ' ' ^ (E^ - 2429. )2 + 0.14823^ ^ 

Here, the energy of photon, E^, is given in unit of keV. In the formulae for the broad positive- 
parity states (1/2+, 5/2+), the energy dependence of the 7-decay width (F^) is explicitly taken 
into account; 

T^ = —a{hc)-'E^B{El), (5) 
9 

where a is the fine-structure constant. In contrast, Eq. (^ results from the energy-independent 
F^ for the sharp 5/2^ state. 
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The total cross section of photodisintegration of ^Be is the sum of three contributions given 

by, 

CT^,n = Cr^,„(l/2+) + 0-^,„(5/2+) + f^SBe (J^,n(5/2"), (6) 

where the prefactor fn»Be is the branching ratio of the 5/2^ state to the n + ^Be decay channel. 



We take a 7 % decay {fn»Be = 0.07) based on the experimental studies |^ |3^ to derive the 
reaction rate. The reaction rate does not change largely even if we assume a 100 % decay 
{fn^Be = 1) to the n + ^Be decay channel. This point will be discussed later in section ^ 

The photoneutron cross section can be converted into the neutron capture cross section of 
^Be based on the detailed-balance theorem. The neutron capture cross section, UnBe, is related 
with the photoneutron cross section, cr^^„, by the formula. 



2(2j, + 1) k^, 
{2jh + l){2j, + l)kl' 



CnBeiEn) = J^. -, wr, ■ T2'^7,ri{Ey) , (7) 



where a, b and c denote ^Be, ^Be and neutron, respectively. The energy, £"„, is the neutron 
energy with respect to the threshold of n+^Be, defined as, 

En = E^- 1.6653, (8) 

in MeV. The wave numbers kn and k^ are given, respectively, by, 

- |. (10) 

where /i is the reduced mass of the n+^Be system. Thus, the neutron capture cross section is 



given by, 

a„Be(K) = 2.415 X io-3 (-^^ + 1-6653)^ ^^^^^^^^^ ^^^^ 

from the measured photoneutron cross section, Eq. (^. Note that (= 10^ E^) in Eq. (|]) is 
calculated as E^ = En + 1.6653, through Eq. (||). 
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4 Theoretical analysis 



4.1 aan reaction rate 

We derive the thermal reaction rate of a{an, 7)^Be in two steps through a metastable ^Be. We 
apply the formulation, which was used to derive the triple alpha reaction by Nomoto et al. [|T^] , 
to the ^Be formation. The same formulation is also employed in the NACRE compilation [|I^]. 
The reaction rate N\ < aan > (i.e. A^^ < av is calculated by, 

Ml < (^) {J;^f f eM-E/k.T)N. < .V >- E,E, 

(12) 

where ^aa is the reduced mass of the a + a system. E is the energy of a + a with respect to the 
threshold (the sum of the rest masses of a + a). aaa{E) is the cross section for ^Be formation 
at an energy E, which is, in general, different from the ground-state resonance energy of ^Be, 
Es^e- ro(^Be, i?) is the a-decay width of ^Be at an energy E. We describe the '^Be formation 
and its a decay in the next subsection. 

The neutron capture rate of ^Be is calculated by, 

Na < av >"'^^= Na (4^) (ttT^Y' r a^s^eiE'; E) exp{-E' /kBT)E'dE', (13) 

where iin^Be is the reduced mass of the n + '^Be system. E' is the energy of n + ^Be with respect 
to the threshold (the sum of the rest masses of a neutron and a ^Be). Note that this threshold 
varies with the formation energy E of ^Be. an»Be{E'; E) is the cross section for neutron capture 
of ^Be at the formation energy E. 



We calculate this cross section using the result of Eq. (|TT|) in section |3.2| . Since the neutron 



energy in Eq. (]TT|) is measured from the threshold for n + Be, we convert E' to En with the 



ground-state resonance energy Es^e by, 

En = E' + E- EsBe, (14) 
to calculate the neutron capture cross section. 
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4.2 ^Be formation 

We calculate the ^Be formation cross section following the prescription of Nomoto et al. For 
details, see Section 2.2.2 of Ref. [|18l. The cross section is given by, 



a^^{E) = ^ exp[-{EG/Ey/% (15) 

where Eg is the Gamow energy given by, 

^ 0.98948^1 ^2^^/^ MeV^/^ (16) 

Here Zi, Z2 are the proton numbers of target and projectile nuclei. A is the reduced mass 
number, 

calculated from the mass numbers of target and projectile nuclei Ai, A2. 
The S-factor S{E) is calculated by, 

^^^^ ^ 0.6566c.. ^^^r^E.)r,^g^^^^ exp[(Ec/i^.)^/^ + - a,E)] MeVbarn, (18) 

where Er is the resonance energy and Ur = [(27. + l)/(2Ji + 1)(2J2 + 1)](1 + ^12)- The energy 
dependent total decay width T{E) is calculated by, 

r{E) = r{Er) expl-iEc/Ey/' - aiE + {Ec/Erf/^ + aiE,]. (19) 

Note that r2(^) in Eq. cancels out with Vc,{^Be,E) in Eq. (|T2p. Taking / = 0, the factor 
ai is defined by 

where Ki{x) is the modified Bessel function of order I, Er is given by, 

2(1 9 

^H = ^MeV (21) 
with R = 1.44(y4j;''^ + A^^^) fm and x is given by 

X = 0.525{AZiZ2RY^^. (22) 
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5 Numerical results 



5.1 Calculated reaction rate 



We numerically integrate Eqs. ([12D and (0) with the cross sections Eqs. (0) and (JTSj) to 
derive the reaction rate A^^ < aan >. We provide the calculated data of the reaction rate 
A^^ < aan > in Table |^. The grid points of temperature is the same as those of the NACRE 
compilation for comparison. 

Figures ^ and ^ display the numerical result of the reaction rate A^^ < aan > as a function 



of temperature (Tg). For comparison, the reaction rates by Caughlan and Fowler (CF88) ||T2 
and Angulo et al. (NACRE) are shown by dashed and dotted curves, respectively. 

The difference from the CF88 is apparent at Tg < 0.028 due to lack of the off-resonant 
contribution in CF88. The neglect of the off- resonant contribution leads to the underestimate 
of the reaction rate below Tg = 0.01 by more than 20 orders of magnitude. The proper treatment 
of the off-resonant contribution from the metastable ^Be nucleus is, therefore, essential at low 
temperatures. The difference remains substantial at high temperatures as seen in Figure ^ 
The CF88 rate is larger than the present rate at 0.028 < Tg < 3 and becomes smaller above 
Tg ~ 3. On the other hand, the NACRE rate is significantly larger than the present rate at 
low temperatures Tg < 0.028 (Fig. ^), whereas it is in good agreement with the present rate at 
high temperatures Tg > 0.1 (Fig. |). 

The ratios of the CF88 rate and the NACRE rate to the present reaction rate are shown 
in Fig. P by the dashed curve and the dotted curve, respectively. The difference of the CF88 
from the present rate amounts to a factor of 2 at Tg > 0.1. This is because the CF88 is based 
on the photoneutron cross section of Berman et al, which shows a larger cross section (1.6 mb) 
peaked at an energy (1.671 MeV) closer to the threshold than the present cross section. At 
Tg > 3, the CF88 rate becomes smaller because it neglects other resonances (mainly the 5/2~ 
state) at higher energies. The NACRE rate is larger by factors 4-12 than the present rate in the 
low temperature regime (Tg < 0.028). It becomes closer to the present rate in the transitional 
regime (Tg ~ 0.028) from off-resonance to on-resonance and consistent with the present rate to 
within ±20% above Tg = 0.1. 
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The discrepancy between the NACRE rate and the present rate at Tg < 0.028 is attributed 
to different behaviors of low-energy tails of broad resonances. Presently, the cross section for 
the 5/2+ state is formulated with the energy-dependent as in Eq. (|^), while the NACRE 
compilation adopted energy-independent widths for the same state (T^ = 0.90 ± 0.45 eV) and, 
in addition, for the l/2~ state (F^ = 1 W.u. = 0.45 eV). As a result, the low-energy tail 
decreases faster with energy in the present cross section than that in the NACRE compilation. 

The 1/2+ resonance immediately above the neutron threshold plays an essential role in the 
reaction rate. The present photoneutron cross section for the 1/2"*" resonance agrees well with 
that of Fujishiro et al. which was employed in the NACRE compilation. This explains a 



similarity between the present rate and the NACRE rate at Tg > 0.1. 

The 5/2~ resonance contributes to the reaction rate only at high temperatures. In the 
current study, we have taken the 7 % decay of the 5/2~ state to the n -|- ^Be channel. To see 
the sensitivity to the branching ratio, we calculated the reaction rate assuming the 100 % decay 
(/nSfie = 1 in Eq. (||)). It was found that the 100 % branching does not change the reaction 
rate below Tg=l and increases the rate only by 5 % around Tg=10. 

5.2 Analytical expressions 

We provide the analytical formula which reproduces the numerical result of the reaction rate 
N\ < aan > given in Table ^ We fit the numerical data by using the same formulae as those 
adopted in Table 3 of the NACRE compilation [|13|. The reaction rate of the ^Be formation, 
Na < aa > formula, is cxprcsscd by 

Na < aa > formula = 2.43 X 10'^T9"'/^exp(-13.490Tg~'/^ - {Tg/0.15)^) x (1 + 74.5 Tg) 

+6.09 X 10^ Tg"^/^ exp(-1.054/Tg). (23) 

The fitting formula of the reaction rate A^^ < aan > for Tg < 0.03 is parameterized as 

Nl < aan >fitting= Na < aa > formula xai(l + aaTg + a^Tg + a^Tg + a^Tg + a^Tg), (24) 

and for Tg > 0.03 as 

Nl < aan >fitting= Na < aa >formuia xai(l + a2 log^o Tg + a3(logio Tg)^ + a4(logio Tg)^), (25) 
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where are coefficients to be determined. We perform a least-squares fit to the numerical data 
of Table The best fit formula for Tg < 0.03 is 

N\ < aan > fuung = Na < aa > formula X 1.376 x 10"^^ 

x(l -58.80 Tg - 1.794 X 10^T2 + 2.969 x 10^ T| 

-1.535 X 10^ Tg"^ + 2.610 x 10^ T|), (26) 

and for Tg > 0.03 is 

Nl < aan > fitting = A^A <«« >/ormnZa x2.630 x 10"^^ x (1 - 1.359 log^oTg 

+1.995 X 10^1 (logioTg)2 + 5.251 x 10"i(logio Tg)^). (27) 

The above formulae reproduce the numerical data of the reaction rate in Table |^ within ~ 20 
% deviation. 



6 Summary and discussions 

We have studied the astrophysical reaction rate for a{an,'-f)^Be, which is important in nucle- 
osynthesis under alpha- and neutron-rich environments. The formation of ^Be by this three 
body reaction plays an essential role to bridge the mass gap at A=8 and to create intermediate- 
to-heavy mass nuclei starting from light elements. This reaction rate influences largely the 
alpha-rich freeze-out followed by the r-process nucleosynthesis in neutrino-driven winds and 
possibly other nucleosynthesis in various astrophysical sites. 

We adopted the latest experimental data on the photodisintegration of ^Be taken with 
laser-electron photon beams. This experiment covered all resonance states of astrophysical 
importance in ^Be. We added data points near the neutron threshold with the flux correction 
for incident photons and updated the least-squares analysis of extracting resonance parameters. 

We numerically calculated the thermal average of cross sections for the formation of ^Be 
via the metastable ®Be, i.e., a + a ^ ®Be(n,7)^Be. We took into account both off-resonant 
and on-resonant contributions from the ground state in ^Be. The off-resonant contribution, 



which was not included in the previous brief report ||T5 on the photodisintegration experiment. 



becomes important at low temperatures below Tg = 0.028. 
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We have provided the reaction rate in the wide temperature range from T9=10~^ to Tg=10^ 
both in a tabular form and in analytic formulae. The calculated reaction rate is compared with 
the reaction rates of the CF88 and the NACRE compilations. The CF88 rate, which is invalid 
at Tg < 0.028 due to lack of the off-resonant contribution, differs from the present reaction rate 
by a factor of 2 at high temperatures Tg > 0.1. The NACRE rate is 4-12 times larger than the 
present reaction rate at low temperature (Tg < 0.028) and becomes closer to the present rate in 
the transitional temperature region from off-resonant to on-resonant processes. The NACRE 
rate is consistent with the present rate at Tg > 0.1 to within ±20%. Because of the systematic 
treatment of photoneutron cross sections obtained from the dedicated experiment, we would 
like to recommend to use the new reaction rate for possible astrophysical problems. 

It is worthwhile mentioning possible impacts of the new reaction rate on modeling nucle- 
osynthesis in various astrophysical sites. The production of seed elements in the neutrino-driven 
wind and the alpha-rich freeze-out scenarios can be enhanced due to the slight increase of the 
reaction rate at Tg > 3.5. This may lead to a slightly smaller neutron-to-seed ratio at the onset 
of r-process, and hence make an appreciable influence on the entropy condition and expansion 
dynamics of the neutrino-driven wind [Q. Kajino et al. have recently studied the sensitivity 
of the r-process yields to the a(an, 7)^Be reaction rate, and found that a factor of two larger 
rate would change the final r-process abundances drastically in rapid expansion models of the 



neutrino-driven wind These possible implications should be studied in detail by using 

the present new reaction rate. It is to be clarified quantitatively if the ^(an, 7)^Be reaction 
path still predominates the seed production and to what extent the final r-process yields may 
change by including the present reaction rate and also newly-identified competing reactions 
7)^Li(n, 7)^Li(a, ?2)-'^^B [^, in the nuclear reaction network. 

It is also postulated theoretically that the production of ^Be in the He layer of AGB stars 
may change the neutron-rich environment required for the s-process nucleosynthesis. Once 
proton mixing into C-rich He layer could occur either in the interpulse phase [^, or in 
the thermal pulse phase fiy] of the AGB stars, the reignition of H-burning proton capture 
on abundant ^^C would form ^^C. Temperature might exceed Tg = 0.09 so that ^^C{a,nY^O 
provides plenty of free neutrons. Since this neutron production could occur in the He layer, the 
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a{an, 7)^Be reaction may follow and play a role of neutron-poison which regulates the s-process 
nucleosynthesis. As the burning time scale of aan depends strongly on the local neutron density 
and temperature, it is highly desirable to calculate precise profile of the neutron exposure and 
the temperature distribution in order to identify above theoretical speculation. Note that, 
should a(an, 7)^Be reaction occur in such an environment of AGB stars, it is hard to observe 
^Be because ^Be is easily destroyed by the {p, a) reaction in the next thermal pulse. 

Another neutron-rich site in cosmological nucleosynthesis is baryon inhomogeneous Big- 



Bang Universe |jTl|. In segregated fluctuations of neutron and proton number-density distri- 
butions it is identified that the ''Li(r;,, 7)^Li(a, ?t,)^^B and ^Li(t, ?T,)^Be reactions are the most 
viable nuclear flows to produce intermediate-to-heavy mass nuclei in the Big-Bang nucleosyn- 
thesis [^, It should be studied theoretically if the a(an, 7)^Be reaction can be an 



alternative flow path or a competing reaction with those identified previously. 

The a(an, 7)^Be stellar reaction rate was evaluated in a more consistent manner, where 
the systematic experimental data of photodisintegration of ^Be was used, than in the preceding 
compilations |[T^, |T^. As of the present evaluation, however, there still remains some uncertainty 
regarding the contribution from the ^He -|- a channel due to a lack of experimental information 
on low-lying states in ^Be. A further experimental investigation with emphasis on the ^He + 
a configuration is advisable. 
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Table 1: Resonance parameters deduced in the present experiment. 



r 


XX 


Er 


B{X\ i) 


r. 


r„ 






[MeV] 


El:[eVm^] 


[eV] 


[keV] 








Ml:[{eh/2McY] 






1/2+ 


El 


1.735 ±0.003 


0.104 ±0.002 


0.568 ±0.011 


225 ± 12 


5/2- 


Ml 


2.43 


0.295 ± 0.072 


0.049 ±0.012 




5/2+ 


El 


3.077 ± 0.009 


0.0406 ± 0.0007 


1.24 ±0.02 


549 ± 12 
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Table 2: The thermal reaction rate N\ < aan > at representative temperatures. 







rate 


Tg 




rate 


0.001 


1 


lOE- 


59 


0.14 


3 


35E- 


08 


0.002 


6 


30E- 


48 


0.15 


4 


96E- 


08 


0.003 


3 


12E- 


42 


0.16 


6 


92E- 


08 


0.004 


1 


21E- 


38 


0.18 


1 


18E- 


07 


0.005 


4 


31E- 


36 


0.2 


1 


78E- 


07 


0.006 


3 


81E- 


34 


0.25 


3 


47E- 


07 


0.007 


1 


36E- 


32 


0.3 


5 


06E- 


07 


0.008 


2 


62E- 


31 


0.35 


6 


29E- 


07 


0.009 


3 


19E- 


30 


0.4 


7 


13E- 


07 


0.01 


2 


75E- 


29 


0.45 


7 


62E- 


07 


0.011 


1 


81E- 


28 


0.5 


7 


83E- 


07 


0.012 


9 


63E- 


28 


0.6 


7 


73E- 


07 


0.013 


4 


30E- 


27 


0.7 


7 


25E- 


07 


0.014 


1 


66E- 


26 


0.8 


6 


62E- 


07 


0.015 


5 


68E- 


26 


0.9 


5 


97E- 


07 


0.016 


1 


75E- 


25 


1 


5 


36E- 


07 


0.018 


1 


30E- 


24 


1.25 


4 


06E- 


07 


0.02 


7 


33E- 


24 


1.5 


3 


12E- 


07 


0.025 


2 


45E- 


22 


1.75 


2 


44E- 


07 


0.03 


2 


98E- 


19 


2 


1 


95E- 


07 


0.04 


1 


36E- 


15 


2.5 


1 


32E- 


07 


0.05 


1 


98E- 


13 


3 


9 


61E- 


08 


0.06 


5 


20E- 


12 


3.5 


7 


39E- 


08 


0.07 


5 


15E- 


11 


4 


5 


94E- 


08 


0.08 


2 


79E- 


10 


5 


4 


20E- 


08 


0.09 


1 


OlE- 


09 


6 


3 


19E- 


08 


0.1 


2 


79E- 


09 


7 


2 


54E- 


08 


0.11 


6 


29E- 


09 


8 


2 


07E- 


08 


0.12 


1 


22E- 


08 


9 


1 


72E- 


08 


0.13 


2 


12E- 


08 


10 


1 


45E- 


08 
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Figure 1: Neutron detection efficiency as a function of neutron energy. 
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Figure 2: Energy distribution of the laser-electron photon beam measured with a high- purity 
Ge detector. The energy distribution, which best reproduces the response of the Ge detector, 
is shown by the solid curve. 
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Figure 3: Photoneutron cross sections for ^Be. A low-energy portion is magnified in an inset. 
The data below 2 MeV are corrected for the 7-ray flux responsible for the photodisintegration. 
The corrected data are shown by the solid circles. The data at the three lowest energies (1.69, 
1.71, 1.73 MeV) near the neutron threshold are consistent with the data of a radioactive isotope 
measurement (slashed squares) |3^ as seen in the inset. The best least-squares fit is shown by 
the solid lines (thick solid line for sum, thin solid lines for breakdown). 



29 




Figure 4: The reaction rate A'^^ < aan > in the present study is shown by the sohd curve as 
a function of temperature. The reaction rates by CF88 [T^ and NACRE ||TB| are shown by 
dashed and dotted curves, respectively. 
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Figure 5: The zoom up of Fig. |^ for a high temperature range. The reaction rates A^^ < aan > 
are shown as functions of temperature between Tg = 10~^ and T9=10 with the same notation 
as in Fig. ^. 
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Figure 6: The ratios of the reaction rate N\ < aan > of the CF88 and NACRE compilations 

to the present reaction rate are shown by dashed and dotted curves, respectively, as functions 
of temperature. 
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